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We use a chiral random matrix model with 2Nf flavors to 
mock up the QCD Dirac spectrum at finite chemical poten- 
tial. We show that the 1/N approximation breaks down in 
the quenched state with spontaneously broken chiral symme- 
try. The breakdown condition is set by the divergence of a 
two-point function, that is shown to follow the general lore of 
macroscopic universality. In this state, the fermionic fluctua- 
tions are not suppressed in the large N limit. 



1. In the presence of a chemical potential, the lattice 
QCD Dirac operator is non-hermitean. As a result, the 
fermionic determinant (probability measure) of the QCD 
Dirac operator is a complex number. Since Monte-Carlo 
simulations demand a positive definite measure, straight- 
forward algorithms have failed Qj. Unconventional QCD 
algorithms have been devised, leading to a chiral phase 
transition at finite chemical potential and strong coupling 
(g = oo) [§]-[|. In the quenched approximation, however, 
the restoration of chiral symmetry was found to set in 
at surprisingly small chemical potentials of the order of 
half the pion mass, thereby vanishing in the chiral limit 
j| . The origin of the discrepancy was traced back to the 
hase of the Dirac operator ||, as demonstrated by in 
using lattice simulations. 
While many of the basic issues related to the quenched 
calculations are known on the lattice ^-||, we feel the 
need for a better understanding using simple models. 
The spontaneous breaking of chiral symmetry in QCD re- 
flects the spectral distribution of quark eigenvalues near 
zero-virtuality [|| . For sufficiently random gauge configu- 
rations, the Dirac operator can be approximated by a chi- 
ral matrix with random entries |^,p| , in qualitative agree- 
ment with a number of lattice simulations for both the 
bulk eigenvalue distributions M , and correlations |l0|,|llj • 
Chiral random matrix models with finite chemical po- 
tential have been discussed recently in [p~2| — p~4[| using var- 
ious methods. For zero chemical potential, the ground 
state breaks spontaneously chiral symmetry. The low- 
lying spectrum is that of constituent quarks with no 
Fermi surface. Naively one would expect that for a chem- 
ical potential of the order of the constituent quark mass, 



chiral symmetry would be restored. This expectation 
is not borne out by direct numerical calculations which 
show that chiral symmetry is restored for any nonzero 
chemical potential in the chiral limit and for large matrix 
sizes [^2) , in qualitative agreement with the quenched lat- 
tice simulations. The numerical results have been shown 
to follow from a modified replica trick |l2|,|l5|, in agree- 
ment with ||. 

In this letter, we reanalyze the spectral density with 
2Nf conventional quarks. In section 2, we introduce 
the model and discuss the 1/N expansion around the 
quenched state with spontaneously broken chiral sym- 
metry. In section 3, we evaluate a pertinent two-point 
correlation in the same state, and show that it diverges in 
the eigenvalue plane in the small eigenvalue region. The 
support for the ensuing spectral distribution is in agree- 
ment with a recent result fH| using different arguments. 
We explain why. Our result for the two-point correlator, 
extends the macroscopic universality argument in ran- 
dom matrix models to the chiral case. In section 4 
we discuss the difference between the unquenched and 
quenched spectral densities. We argue that the quenched 
state with spontaneous symmetry breaking does not sup- 
press the fermionic fluctuations. Our conclusions are 
summarized in section 5. 

2. A model form of the QCD Hamiltonian in the chi- 
ral basis with finite chemical potential fi and zero current 
quark mass is 
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where only constant quark modes have been retained [ p"3"[ . 
In the context of four dimensional field-theories Q5 is 
the Hamiltonian of a Dirac particle in five dimensions. 
The translation of this operator to chiral random matrix 
models is achieved by assuming that the gluonic configu- 
rations A are sufficiently random so that A — > R, where 
R is an N x N complex random matrix with Gaussian 
weight. With this in mind, (Q) is the sum of a determin- 
istic D and random (chiral) piece R. The deterministic 
piece is non-hermitian for /1 / 0, with complex ±i|/i| 
eigenvalues. 

Now, consider the partition function associated with 
(0) for 2N f flavors 
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Z(2N f ,fi) =< det 2N f(z - M) > 



(2) 



with V2N f = —log Z/2Nf playing the role of a complex 
potential. For z = im, (0) mocks up the partition func- 
tion of 2Nf quarks of equal mass m. The averaging in 
(0) is carried using a gaussian weight 
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with E setting the scale of chiral symmetry breaking 
(throughout S is set to 1 for simplicity, unless indicated 
otherwise) . The gaussian weight (||) is an idealization of 
the gluon action, and defines the quenched averagings to 
be discussed below. For sufficiently small quark eigenval- 
ues (here constant modes), the gluon interaction may be 
thought as random |?]] . The gaussian weight follows from 
the principle of maximum entropy. 

First, let us consider the spectral distribution associ- 
ated to (H|) in the quenched approximation. With V2N f 
playing the role of a complex potential, the spectral dis- 
tribution follows from Gauss law [fl5[ 

v(z, z) = lim p(z,z)= lira -—r:d z d z V2Nf (4) 

Explicitly, ■nv{z 1 z) = dG/dz, with the resolvent 



G(z,z) = ±(Tr(z-M)- 1 ) 



(5) 



We have retained (z, z) in both (|j) and @ to allow for 
the possibility that holomorphic symmetry may be bro- 
ken in the thermodynamical limit by the quenched aver- 
aging (|J). The factor of 1/N in (j|) implies that v(z,z) 
is normalized to one on its pertinent support. 

For large z and to leading order in 1/N, G obeys Pas- 
tur's equation (planar approximation) |l3|.[L9[ , 



G((z-G) 2 +A* 2 )~z + G = 



(6) 



The solution to (g|) is a holomorphic function of z, ex- 
cept for line discontinuities, with end-points given by the 
zeroes of the discriminant 

4/iV + z 2 (8^ 4 - 20^ 2 - 1) + 4(^ 2 + l) 3 = (7) 

For [i — 0, there is a cut along the real axis be- 
tween —2 and +2. The discontinuity along the cut is 
Wigner's semicircle, with a quenched and stable state 
that breaks spontaneously chiral symmetry. With in- 
creasing p 2 — > 1/8, the two roots approach each other, 
followed by the emergence of two new roots from real in- 
finity. For /i 2 = 1/8 they coalesce pairwise, and move to 
the complex plane for p 2 > 1/8. This behavior suggests 
that the quenched state with v(+i0, —iQ) ^ 0, supports 
a chiral condensate up to p 2 = 0.125, after which a first 
order transition is observed. This conclusion, however, 
does not hold since the 1 /N approximation breaks down 
for small eigenvalues in the present model, as we now 
show. 



3. The breakdown can be probed either by looking 
at the nonplanar corrections to (^J), or by testing the 
the stability content of the quenched state against local 
correlations. One pertinent correlation function is given 
by the two-point function 



iV 2 C(z, z) =< Tr(z - M) _1 Tr(^ - M^y 1 > c 



(8) 



in the quenched state. The lower script in (|J) is for con- 
nected. The physical interpretation of (|J) will be clarified 
in the next section. For large values of z, (pi) may be an- 
alyzed using again a 1/N expansion [ Lq-jl8[. Following 
the arguments in |l6|], the result is [^1 , 



N 2 C(z,z) 



^3 2 Log{[(H- M 2 ) 2 



G| 4 ]/H 2 } (9) 



where H — \z — G| 2 /|G| 2 and G is a solution to Pas- 
tur's equation (|6[). For fj, = and z = z outside the 
cut, (||) reduces to N 2 C{z 7 z) = l/z 2 (z 2 - 4) 2 , which is 
different from the hermitean (nonchiral) result of [p.6| , 
but in agreement with the one of p2] |. The result (Q) 
shows that in our case the two-point correlation function 
is also amenable to the one-point function, thereby ex- 
tending the macroscopic universality argument discussed 
originally in jl6) to the chiral case. 

From (j^) we observe that the two-point correlation 
function diverges in the eigenvalue-plane in the domain 
prescribed by the zero of the logarithm, 



G| 2 (1-|G| 2 )- M 2 |G| 







(10) 



For /i = 0, this condition is fulfilled for G = z, which is 
a trivial (unphysical) solution to Pastur's equation, and 
for |G| 2 = 1 a nontrivial (physical) solution to Pastur's 
equation along the cut. The latter when supplemented 
with the condition lirri^oo G ~ 1/z in the outside, yields 
Wigner's distribution for the eigenvalue density. 

Figure |l| shows the envelope in dashed line for which 
the condition (^) is met in the w = iz plane, for several 
values of the chemical potential /i. A structural change 
occurs at /i 2 = 1. The solid lines are the mean- field 
solution following from (Q). They lie in the eigenvalue- 
domain where the 1/N fluctuations are dominant, a sig- 
nal that chiral symmetry is restored. Since (||) mixes z 
and z, holomorphic symmetry is spontaneously broken in 
the quenched state. One can recast (|l0|) using (j|) into 
(iw = x + iy) 

( M 2 -x 2 )V = [4 M 4 (l-^ 2 )-(l^ 2 -8 M 4 )x 2 -4 At 2 a ; 4 ] , (11) 



The result Qll| ) is in agreement with the one in JL2| , where 
the quenched spectral density was evaluated in the inside 
(small eigenvalues) using a pair of conjugate quarks in 
the quenched approximation (essentially \z — M| in (g)). 
Conventional quenched QCD in the 1/N approximation 
works as well outside. The two approaches agree on ( |Tl| ) 
which sets their domain of validity. To be able to de- 
scribe the spectral distribution inside using the quenched 
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version of (0) requires a different method than the l/N 
approximation we have used. 







FIG. 1. The envelope (dashed lines) is from eq. ( Jlo| ) in the 
plane w = iz while the cuts (solid lines) are from eq. (Q), for 
different y. Shaded islands represent "other vacuum". 

Figure ^ shows the analytical behavior (solid lines) for 
the imaginary part of the resolvent (§]§) and the corre- 
lation function (|^) for two values of w = iz and fixed 
/i 2 = 2, that is w = iA — y (left) and w — i0.02 — y 
(right) with y > 0. The dashed curves are a comparison 
to a numerical estimate using an ensemble of 200 chi- 
ral gaussian plus deterministic matrices with dimension 
100 x 100. Note that in the right figures, two different 
solutions to (^J) have to be used while crossing the re- 
gion of non-analyticity (shaded region). In the "bounded 
region" between the islands of Figure [l] the asymptotic 
condition lim z _ +00 G ~ 1/z is no longer required. We 
remark that ImG = for z = 0, and chiral symmetry is 
restored. 
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FIG. 2. Imaginary part of the resolvent (upper figures) and 
the correlation function (lower figures) for w = ix — y with 
x = 4, 0.02 and y > at y 2 = 2. The solid curves are 
analytical while the dashed curves are numerical (see text). 

4. A dequenching of the spectral density (|J) can be 
simply achieved by noting that since (z — M) is non- 
hcrmitean, we can split it into a phase and a modulus 
through 



(z - M) 2 = \z- M| 2 x (- 
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In terms of (^|) and (IT4), the unquenched spectral density 
for finite N and Nf is now given by p(z, z). From (||) the 
contribution of the modulus is 
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Pmod(z, z) = V(Z, Z) 7^C(Z, z) 
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(13) 



while the contribution of the phase is 
NfN 



P P h(z,z) 
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The density p(z, z) is the sum of |[|) and (|l|) modulo 
an extra contribution from the crossed terms, essentially 
the last term in (|l^) to the order quoted. We note that 
the connected two-point function (^) is essentially the 
fluctuation of the phase of the fermion determinant in 
the quenched state. It appears with opposite signs in the 
modulus and the phase, and cancels in the sum. The dif- 
ference p{z, z) — v(z, z) accounts for the effects of the sea 
fermions (unquenching) . This difference involves nonlo- 
cal correlation functions of z and z, of which (Q) is a 
generic example. The fermionic effects are down by 1/./V, 
provided that the correlation functions are stable. This 
is not the case in the quenched state with spontaneous 
symmetry breaking. 

Further insights to the above results can be achieved, 
if were to note that (|J) may be rewritten as 



N 2 C(z,z) = (qUQ^Q)c 



with 
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(15) 



(16) 



playing the role of "propagators". For [i = 0, (q^q) and 
(Q^Q) are nonzero for z ~ iO, and the vacuum contribu- 
tions cancel out. Typically, p(z, z) — is(z, z) ~ 0(Nf /N) 
for z ~ iO. For fi = and large N, the fermionic fluc- 
tuations are suppressed in the quenched state. The lat- 
ter breaks spontaneously chiral symmetry, and preserves 
holomorphic symmetry (q «-> Q). 

For small /1, ( |l5| ) through ( |lC)| ) diverges when closing 
on the dashed curve of Figure 1 from the outside. This 
is a signal that ( |l5| ) is receiving increasingly large con- 
tributions from the "mixed condensates" (q^Q) and their 
conjugates (other "vacuum"), so the fluctuations are not 
suppressed any more. As a result p(z, z) — v(z, z) ~ 
O(NfN ) for small z. In the quenched state with sponta- 
neous symmetry breaking, the fermionic fluctuations are 
not suppressed in the large N limit. 

5. Using a chiral random matrix model with a fi- 
nite chemical potential /1, we have shown that the small 
and quenched eigenvalue distribution of the Dirac op- 
erator is fluctuation driven, with a size in the complex 
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plane conditioned by the divergence of the connected 
two-point correlation function (|8j). In this domain the 
X/N approximation breaks down and the effects from 
the fermionic fluctuations do not decouple in the ther- 
modynamical limit. The fluctuation driven phase breaks 
spontaneously holomorphic symmetry and preserves chi- 
ral symmetry, as originally shown in Ji~2j |. 
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